Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 2 February 2008 



(MN M?eX style file vl.4) 



Spectral lags and the energy dependence of pulse width in 
gamma-ray bursts: contributions from the relativistic 
curvature effect 

Rong-feng Shen Li-ming Song, and Zhuo Li 

Particle Astrophysics Center, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100039, China 
Accepted 2005 April 27. Submitted 2004 July 20 

ABSTRACT 

We compute the temporal profiles of the gamma-ray burst pulse in the four Burst and 
Transient Source Experiment (BATSE) Large Area Detector (LAD) discriminator 
energy channels, with the relativistic curvature effect of a expanding fireball being 
explicitly investigated. Assuming an intrinsic "Band" shape spectrum and an intrinsic 
energy-independent emission profile, we show that merely the curvature effect can 
produce detectable spectral lags if the intrinsic pulse profile has a gradually decaying 
phase. We examine the spectral lag's dependences on some physical parameters, such 
as the Lorentz factor F, the low-energy spectral index, a, of the intrinsic spectrum, 
the duration of the intrinsic radiation t'^ and the fireball radius R. It is shown that 
approximately the lag oc and oc t^, and a spectrum with a more extruded shape (a 
larger a) causes a larger lag. We find no dependence of the lag on R. Quantitatively, 
the lags produced from the curvature effect are marginally close to the observed ones, 
while larger lags require extreme physical parameter values, e.g., F < 50, or a > —0.5. 
The curvature effect causes an energy-dependent pulse width distribution but the 
energy dependence of the pulse width we obtained is much weaker than the observed 
W (X E~^-^ one. This indicates that some intrinsic mechanism(s), other than the 
curvature effect, dominates the pulse narrowing of GRBs. 
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; 1 INTRODUCTION 

, Cheng et al. (1995) were the first to analyse the spectral 
lag of gamma-ray bursts (GRBs), which is the time delay 
between the peaks in the Burst and Transient Source Ex- 
periment (BATSE) Large Area Detector (LAD) Channel 1 
(25 - 50 keV) and Channel 3(100 - 300 keV) light curves. 
Subsequently, several authors have carried out more analy- 
sis work on the GRB lags. Norris et al. (1996) and Norris, 
Marani & Bonnell (2000) found the cross-correlation func- 
tion lags between BATSE Channel 1 and Channel 3 pho- 
tons tend to concentrate near < 100 ms; for six bursts with 
known redshift, the z-corrected lags are distributed between 
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6 and 200 ms. Wu & Fenimore (2000) extended the analysis 
to very low energy (~ 2 keV); they found that about 20% 
of GRBs have detectable lags and that GRBs don't show 
larger lags at lower energy. A recent measurement by Chen 
et al. (2005) for the BATSE bursts shows that the majority 
of the lags are below ~ 200 ms, and that the histogram of 
the lags peaks around 30 ms. More intriguingly, Norris et 
al. (2000) found that, for those six bursts with known z, the 
peak luminosity is anti-correlated with the lags. This rela- 
tionship provides a useful tool to estimate the distances of 
large sample of GRBs by analysing their light curves. 

Three theoretical explanations for the lag/luminosity 
relation have been proposed: the relationship is due to the 
variation in line-of-sight velocity among bursts (Salmonson 
2000); it is caused by the variation of the off-axis angle when 
viewing a narrow jet (loka & Nakamura 2001); or it is caused 
by radiation cooling - highly luminous burst cools fast and 
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the lag will be short (Schaefer 2004). However, the problem 
of what mechanism(s) causes the spectral lags of GRBs re- 
mains umesolved. Salmonson (2000) did not explain the ori- 
gin of the lag but assumed that it derives from some proper 
decay time scale At' in the rest frame of the emitter. In 
the model by loka & Nakamura (2001), the lag is caused by 
the far side of the emitting region producing lower-energy 
radiation after a longer light-travel time, for a narrow jet 
with viewing angles outside the cone of jet. However, their 
model requirements seem too stringent (sec below). Schae- 
fer (2004) proposed radiative cooling as the origin of the lag. 
The difficulty of this explanation is that in order to adjust 
the observed synchrotron cooling time scale to be compara- 
ble with the lag time scale, the magnetic field has to be ~ 7 
Gauss, a value much below the strength required by most of 
the current models (e.g. Piran 1999). 

Kocevski & Liang (2003) have assumed that the ob- 
served lag is the direct result of spectral evolution, another 
property of GRBs (Norris et al. 1986; Bhat et al. 1994). In 
particular, as the peak energy of the GRB's vF^, spectrum 
decays through the four BATSE channels, the photon-flux 
peak in each individual channel will be shifted, probably 
producing the measured lag. From a sample of 19 GRBs, 
Kocevski & Liang (2003) found an empirical relation be- 
tween the peak energy's decaying rate and the GRB lag. 

It is widely accepted that the gammarrays come from 
a relativistically expanding fireball surface with Lorentz 
factor r > 100 (Lithwick & Sari (2001) and references 
therein). At some distance from the central source (e.g. 
R = 10^2 ~ 10"cm, cf. Piran (1999)), photons emitted 
from the region on the line of sight and those from the side 
region at an angle of S ~ l/F with respect to the line of sight 
are Doppler-boosted by different factors and travel difi'erent 
distances to the observer. This is what we call the curva- 
ture effect. Comparing the radiation from the side region 
and that from the line-of-sight region, for the latter its pho- 
tons are Doppler-boosted to higher energies and arrive to 
the observer earlier; its observed temporal structure will be 
boosted to be narrower. 

The motivation of this paper is to see how the curva- 
ture effect will change the intrinsic pulse profile in different 
energy channels, and special interest is focused on whether 
merely the curvature effect can produce the spectral lags of 
the pulses. Except for the soft photon lags, the pulses in 
GRBs show another temporal property, i.e. the pulse nar- 
rowing at higher energy, or pulse width as a function of 
energy (Fenimore et al. 1995, Norris et al. 1996). We are 
also interested in probing the contributions of the curvature 
effect to these properties. 

A fireball internal-external shocks model has emerged 
for the theoretical understanding of the origin of GRBs (Pi- 
ran 1999, 2004). According to this model, GRBs are pro- 
duced when an ultra-relativistic outflow dissipates its kinetic 
energy through the internal collisions within the outflow it- 
self. The afterglow occurs when the flow is decelerated by 
shocks with the circumburst medium. This model has made 
many successful explanations and predictions to the obser- 



vations of GRBs. Our investigation to the curvature effect 
will be based on the frame of the internal shock model, where 
the shock is generated from the colliding shells, and electrons 
are accelerated in the shock and radiate. 

loka & Nakamura (2001) proposed a model in which 
a narrow jet is viewed at off-axis angle to explain the 
lag/luminosity relation and the variability/luminosity relar 
tion. Their model successfully reproduces the lag/luminosity 
relation, while the lag is caused by the curvature effect of 
the jet, which increases with the off-axis angle. However, 
one substantial problem with this model, £is pointed out by 
Schaefer (2004), is that it works with an exacting assump- 
tion that the jet opening angle always equals F~^. Further- 
more, loka & Nakamura (2001) only consider an instanta- 
neous emission in the jet rest frame, which is a too sim- 
ple assumption. Consequently, in their model there is no 
lag when the jet is viewed near-on-axis. Different from loka 
& Nakamura (2001), we have considered various rest-frame 
emission profiles and assumed an isotropically expanding ra- 
diation surface. 

We present our model in Section 2, including the basic 
assumptions and formulas. The major results are presented 
in Section 3, based on which we give our conclusions and 
discussion in Section 4. 



2 MODEL 

2.1 Three time scales 

Three time scales are involved in determining the tem- 
poral structures of pulses in GRBs: (i)cooling time scale; 
(ii)hydrodynamic time scale; (iii)angular spreading time 
scale (Kobayashi et al. 1997, Wu & Fenimore 2000). 

In the synchrotron cooling model, the shock-accelerated 
electrons cool via synchrotron emission, and the electron's 
average energy becomes smaller and the radiated power de- 
cays. As pointed out by Wu & Fenimore (2000) , the stan- 
dard internal-shock model gives an observed synchrotron 
cooling time scale at a given photon energy as 



huob 



IQQkeV ' 



-1/2 



where es is the equipartition parameter for the ratio of the 
magnetic energy density to the total internal energy density; 
for a typical value in the internal-shock model, e.g. es = 0.01 
{B ~ 10* G), the cooling time scale is far shorter than the 
lag time scale. 

The hydrodynamic time scale is related to the energiz- 
ing of the electrons. In the internal-shock model, if one as- 
sume the the local microscopic acceleration of electrons is in- 
stantaneous, then the hydrodynamic time scale is attributed 
to the shell-crossing time of the shock, T'^y„ = A'/«^^, where 
A' is the shell width and is the shock velocity, both in 
the comoving frame of the upstream flow. We hardly know 
about A'. However if one assume the shells axe radially ex- 



© 0000 RAS, MNRAS 000, 000-000 



Spectral lags in gamma-ray bursts 3 



panding, in the observer frame this time scale is 



' 1015cm ^nOO^ 

where R is the radius at which the shell radiates, and F 
is the Lorontz factor of the shell (Rydc & Pctrosian 2002). 
(Note that apart from the shock acceleration scenario, there 
can be other particle energizing mechanisms, e.g. magnetic 
field rcconncction (Stern 1999), which may have a different 
hydrodynamic time scale) 

The angular spreading time scale is the delay between 
the arrival times of the photons emitted at the line-of-sight 
region and of that emitted at the side region of the shell (Sari 
& Piran 1997). Because of the rolativistic beaming of the 
moving radiating particles, only the emission from a narrow 
cone with an opening angle of ~ l/F is observed. This gives 
a time scale of the delay of 



= 1.7s ( 



^ )(— )-^ 



2.2 Assumptions 

In this paper, different from Schaefer (2004), we assume 
the cooling time scale is much shorter than the other two 
time scales, i.e. the accelerated particles radiate their en- 
ergy rapidly. Thus the rest frame duration of the emission 
is determined by the hydrodynamic time scale. 

We consider a thin shell expanding with a relativistic 
speed, whose Lorentz factor is F. The shell begins to radi- 
ate at radius R. In the co-moving frame of the shell, the 
radiation intensity of the shell surface I^v' ,t') is assumed 
to be isotropic, and has an energy-independent time history 
f'{t'). Note that all the quantities in the rest frame of the 
radiation surface arc labeled with a prime note. 

Band et al. (1993) found that the GRB spectra are well 
described at low energy by a power law with an exponential 
cutoff and by a steeper power law at high energy. The typical 
fitted value distributions for the low-energy spectral index 
(a) is -1.5 ~ -0.3, the high-energy spectral index (/3) is - 
2 ~ -3 and the peak energy (Ep) of the vF,^ spectrum is 
100 ~ 500 keV (see Preece et al. 2000). By modeUng, Qin 
(2002) showed that the relativistic expanding of the fireball 
does not alter the shape of GRB's rest frame spectrum, but 
only shifts the peak of the spectrum to a higher energy; the 
distribution of Ep is scaled with the Lorentz factor F, as 
shown in their Table 4. So we directly adopt the spectral 
form used by Baud ct al. (1993) as the rest frame spectrum, 
and choose values for the low-energy spectral index a to be 
-0.8 and the high-energy spectral index /? be -2.4; The peak 
energy of the rest frame spectrum E'p is adjusted such that 
Ep is 350 keV, where the relation Ep = 1.65 x TE'p, derived 
from Table 4 in Qin (2002) in the range of 10 < F < 2000, 
is used. 



2.3 The formulas 

The radiation intensity in the observer's frame is /(/Lt, v, t) = 
I{n,v)f{t). It is connected with the rest frame radiation 
intensity by 

I{l,,v)f{t) = V\y)l'{u')f{Vt') (1) 
and 

u = v'V{^), (2) 

where fi = cosO, and is the angle of the concerned local 
radiation surface to the line of sight; ©(/i) = [r(l — /i/3)]~i is 
the local relativistic Doppler factor respect to the observer, 
where F is the Lorentz factor of the radiation surface and 
/3 = Vl-F-2. 

We use t to refer to the photon emitting time and use T 
to refer to the time that the photon arrives at the observer. 
We define that, at the time t = 0, the first photons are 
emitted from the surface and, for simplicity, we also tune 
the arrival time T of the first photon emitted from the line- 
of-sight region {jj, = 1) to be too. Then the arrival time of 
a photon emitted from region n at time t is 

T={l^l3)t + {1- y){R + I3ct)/c. 

The first part of the right-hand side of the equation is caused 
by the motion of the shell, and the second part is due to the 

difference between the light travel distances of photons omit- 
ted from the line-of-sight region and from the side region. It 
can be rewritten as 

T=(l-M/3)*+(l-M)r, (3) 

where r = R/c, which connects the photon's arrival time, 
emitting time and emitting place in one equation. 

At time T, the observed flux comes from the photons 
emitted at the region 1 > ^ > //(T, t = 0), where the bound- 
ary n{T,t = 0) = 1 — T/t, calculated from equation (3), is 
the place whose first photon is emitted with arrival time T. 
The observed specific flux can be obtained by integrating 
the radiation intensity over this region 

r-l-T/T 

FAT) = - 



h{fi,t)fiR^{t)diJ, 



(4) 



For a BATSE LAD discriminator channel {va,Vb), the 
observed photon counts flux is 



nab{T) = I 



-dv. 



From equation (4), substituting the emitting time t in the 
integral with T and /i using equation (3), nab{T) can be 
rewritten as 

nab{T) = -[r + (T-T)l3f J 

where Rit) = c[r + f(T, [i)S\ and equations (1), (2) are used. 
So we get the explicit formula that is used in our calculation. 
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2.4 Intrinsic time profile of the radiation 

Here we introduce throe cases of the intrinsic emission time 
profile. 

(i) Rectangulax profile. During a finite duration, the 
emissivity is constant, with the instantaneous rising phase 
and decaying phase: 



m 



1, 0<t<td 
0, t>td 



In our assumptions td is related to the hydrodynamic time 
scale discussed above, but note that ta is defined in the 
observer's frame, without taking into account the Doppler 
boosting (i.e., it is not the observed time scale). In addition, 
we refer to t'd as the td measured in the comoving frame of 
the shell. 

(ii) One-sided exponential profile. The emissivity decays 
exponentially after an instantaneous rising phase: 



0, 

exp[ 



-{t/td)\, 



t < 
i > 



(iii) Symmetric Gaussian profile. The emissivity has 
both a finite-time rising phase and a finite-time decaying 

phase: 



m = 



0, 

exp[-(*: 



t < 

t>0 



where we introduce the coefficient -1.5 in order that the 
emission starts at t = with the radiation intensity of a 
tenth (e"^'^ ^ 0.1) of its peak value. 



3 RESULTS 

3.1 Observed temporal profiles in different energy 
bands and time lags 

Using the equations and the typical parameters introduced 
above, we calculate the observed pulse profiles in the four 
BATSE LAD energy channels, assuming three different in- 
trinsic pulse profiles. The results are plotted in Figure 1, 2 
and 3, respectively. 

First, we calculate the pulse light curves in the four en- 
ergy channels for the rectangular intrinsic radiation profile. 
It has a steady rising phase, followed by a distinct peak, 
as shown in Figure 1. The rising of the flux is due to the 
expanding of the radiation surface. The rising is steady be- 
cause the radiation intensity is constant during a finite time. 
The decay phase of the observed pulse is due to the angular 
spreading effect. The peak occurs when the intrinsic radia- 
tion begins to cease. For various parameter spaces (duration 
of the emission td : 10^ - 10^ s, radius R : lO" - 10^^ cm, 
Lorentz factor F : 50 — 500), we do not detect any time lag 
between the peaks observed at different energies in this case. 

As for the one-sided exponential and the symmetric 
Gaussian emission profiles, their observed light curves cal- 
culated in the four energy channels are illustrated in Figure 




Observer's time 

Figure 1. Observed pulse from intrinsic rectangular time profile 
in four BATSE LAD energy channels. Solid line: 25 - 50 keV; 
dashed line: 50 - 100 keV; dot-dashed line: 100 - 300 keV; dotted 
line: 300 - 1800 keV. The vertical coordinate is in arbitrary units. 
The peaks of the pulses at 4 energies arrived simultaneously. F = 
100, a = -0.8, /3 =- 2.4, R = 3 X lO" cm, td = 0.9 x lO" s. The 
rising phase comes from the expanding of the radiation surface 
and the rising time is determined by duration of the intrinsic 
radiation. The decay phase is due to the angular spreading effect. 



2 and 3, respectively. For these two emission profiles, their 
peaks are more gradual. The maximum of the photon flux 
at higher energy arrives early than at lower energy. If we de- 
fine the time lag as being the difference between the arrival 
times of the peaks in energy Ghannel 1 and 3 or 4, the re- 
produced lags are quantitatively comparable to the observed 
ones (e.g., ~ 10"^ s if corrected for the cosmological time 
dilation, cf. Norris et al. (2000)). Note that here we do not 
take into account the redshifts, z, of the GRBs, which must 
make the observed lags (1 -|- z) times larger. 

The intrinsic rectangular profile cannot produce the ob- 
served lags because in this case the radiation diminishes 
immediately. In the observed pulse, the rising phase comes 
from the expanding of the radiation surface and the ris- 
ing time is determined by duration of the intrinsic radiation 
(~ ^d/(2r^); also see Qin et al. 2004). The decay phase is due 
to the angular spreading effect. Hence if the intrinsic radia- 
tion switches off immediately, the transition from the rising 
phase to the decay occurs abruptly and induces a sharp peak 
in the observed pulse. For this case, the immediate switch- 
off of the intrinsic radiation dilutes the relativistic curvature 
effect, and hence does not produce the peak lags. 



3.2 The pulse width - energy relation 

For the three intrinsic emission profiles, we find that the 
pulse observed in high-energy channel is narrower than in 
lower-energy channel, which is manifested by the fact that 
the FWHMs of the pulses in separated energy channels de- 
crease with the energy in a power-law form. However, the 
pulse narrowing we obtained is less prominent than that ob- 
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Figure 2. Observed pulse from one-sided exponential intrinsic 
time profile in four BATSE LAD energy channels. Solid line: 25 — 
50 keV; dashed line: 50 - lOOkeV; dot-dashed line: 100 - 300 keV; 
dotted line: 300 - 1800 keV; triple- dots- dashed line: the intrinsic 
emission profile. The vertical coordinate is in arbitrary units. T = 
400, a = -0.8, P = -2.4, R = 5x lO" cm, ta = 5 X 10* s. Lagis 
0.03 s. 
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Figure 4. Symmetric Gaussian intrinsic pulse: lag dependence 
on the Lorentz factor, a = -0.8, /3 = -2.4, R = 5 x 10^^ cm, t^ = 
20 s. Power-law fit gives index -0.87 and -0.93, respectively. 

by ^ 6%, hcricc hardly changes the slope of the pulse width 
versus the energy. 
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Figure 3. Observed pulse from symmetric Gaussian intrinsic 
time profile in four BATSE LAD energy channels. Solid line: 25 — 

50 kcV; dashed line: 50 - 100 kcV; dot-dashed line: 100 - 300 keV; 
dotted line: 300 - 1800 kcV; triple- dots- dashed line: the intrinsic 
emission profile. The vertical coordinate is in arbitrary units. F = 
400, a = -0.8, /? = -2.4, _R = 5 x lO" cm, td = 2x 10* s. Lagia =i 
0.005 s. 

served in real GRBs for which pulse width decay power-law 
index ~-0.4; while the pulse width decay index we obtained, 

for instance in the case of Figure 2 (one-sided exponential 
decay emission profile), is -0.13. 

Other than the "Band" spectrum, we also used an al- 
ternative function - a low-cnorgy power law plus the high- 
energy exponential cut-off at E'p - for the rest-frame spec- 
trum. For the typical values of the parameters we have used, 
this changing of spectrum only narrows the Channel 4 pulse 



3.3 Lag's dependence on other physical 

parameters 

The spectral lag is an important observational property of 
the pulse in GRBs in that it may be used to derive the cos- 
mological distribution of GRBs (Norris 2002) and to discrim- 
inate the internal shock signature and the external shock 
signature in the pulses (e.g., Hakkila & Giblin 2004). This 
motivates us to probe the dependences of the peak lags on 
other physical parameters of the simple model. We choose 
the symmetric Gaussian profile as the intrinsic omission pro- 
file, which includes an intrinsic rising phase. We alter the 
Lorentz factor F, the spectral parameters of the rest frame 
emission {a, f3 and E'p), the radius of the radiation surface 
R, and the rest-frame duration t'^ of the intrinsic radiation, 
respectively, and sec how the Channel 1/3 and the Channel 
1/4 peak lags vary with these changes. 

3. 3. 1 Lorentz factor 

Figure 4 shows that the lag decreases with the Lorentz factor 

following Lag oc F"^. We think this is a natural outcome of 
the relativistic boosting of the time structure. Those pulses 
whose lags are larger may come from colliding shells with low 
Lorentz factors, according to the current standard models 
(e.g., Piran 1999). 

3.3.2 Spectral parameters 

Figure 5 shows that the lag increases with the low energy 
spectral index a of the rest frame spectrum of the pulse. A 
larger a means more photons arc concentrated around the 
peaJi energy E'p of the v' Fl spectrum. Then, for a narrower 
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-0.5 
Low energy spectral index a 

Figure 5. Symmetric Gaussian intrinsic pulse: lag dependence 
on a, the low-energy spectral index of the adopted spectrum. F 
= 100, e'j,= 1.75 keV, ^ = -2.4, R = 5x lO" cm, t'^ = 40 s. Note 
that GRB spectrum with a > is very rarely observed (Preece et 
al. 2000); the calculated data points there are only used to show 
the tendency. 
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Figure 6. Symmetric Gaussian intrinsic pulse: lag dependence 
on the observed spectral break energy Ep of the pulse spectrum. 
r= 400, a= -0.8, 13= -2.4, R= 5x10"* cm, t'^= 100 s. Note that Bp 
is obtained through a simple scaling relation, Ep = 1.65 X Fi/p, 
between Ep and the break energy, Ep, of the rest- frame emission 
spectrum, as a result of the simulation by Qin (2002) . 



spectrum, the curvature efltect will work more effectively in 
producing the spectral lags. This conjecture is supported 
when we alter the high-energy spectral index /?. We found 
that a steeper high-energy power-law spectrum produces a 
larger lag, e.g., the Channel 1/4 lag has a 14% increase for 
f3 changing from -2.4 to -3.0. Compared with the Channel 
1/4 lag, the Channel 1/3 lag has a weaker dependence upon 
/3. 

The lag's dependence on the observed break energy of 
the spectrum Ep is shown in Figure 6. The lag has its max- 
imum when Ep falls near the starting energy of the corre- 
sponding high-energy channel that is used in measuring the 



lag (i.e., ~ 100 keV for Channel 3, ~ 300 keV for Channel 
4). 

The above findings about the lag's dependences on the 
spectral parameters are qualitatively consistent with the 
tendency observed in those long-lag wider-pulse bursts by 
Norris et al. (2005). They found that their long-lag (mea- 
sured for Channel 1/3) burst sample has, on average, lower 
Ep (centered around ~ 110 KeV), larger a (harder low- 
energy power law) and smaller (3 (softer high-energy power 
law) , than the bright burst sample analyzed by Preece et al. 
(2000). 

Substituting the "Band" spectrum with an alternative 

one of a single power law plus an exponential high-energy 
cut-off causes no changes to the Channel 1/3 lag, while the 
Channel 1 /4 lag has a ~40% increase if the observed cut-off 
energy Ep is below 300 keV; for Ep > 300 keV, the increase 
of Channel 1/4 lag is much smaller. 



3.3.3 Duration of the emission 

We find that longer rest-frame duration (t'j) of emission will 
cause larger lags, as is shown in Figure 7. This result may 
be associated with an observed tendency that wider pulses 
exhibit longer lags (Norris et al. 1996; Norris, Scargle & 
Bonnell 2001; Norris et al. 2005). 



3.3.4 Radius of the radiation surface 

The lag appears to be independent on R, the radius of the 
radiation surface (see Figure 8). We know R determines the 
angular spreading time scale, Tang ~ -R/(2r^c). This re- 
sult suggests that though in our model the angular spread- 
ing effect is a necessity in causing the lags, the peak lag in 
the pulses is not correlated with the angular spreading time 
scales. 

In addition, we find that the decrease of pulse width 
with the photon energy is dependent upon the low-energy 
spectral index of the radiation spectrum, as shown in Figure 
9. 



4 CONCLUSIONS AND DISCUSSIONS 

By assuming an intrinsic "Band" -shape spectrum and an ex- 
ponential or Gaussian emission profile, we show that merely 
the curvature effect produces detectable soft lags in the GRB 
pulses. Therefore the soft time lags can be a signature of the 
relativistic motion occurring in GRBs. 

The observed Channel 1/3 lags are typically distributed 
among 10~^ — lO"'^ s (Norris et al. 2000). For typical physical 
parameters, i.e., F w 100, t'^ ~ 40 s, a ~ -1, as shown in 
Figure 4-8, the lags produced by the relativistic curvature 
effect are slightly above 10~^ s, marginally close to those 
observed, after considering the cosmological time dilation 
effect if a GRB redshift of 2 is assumed. To account for those 
observed larger lags (~ 0.1 s), it requires extreme physical 
parameter values, e.g., F < 50 or a > -0.5. 
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Figure 7. Symmetric Gaussian intrinsic pulse: lag dependence 
on t^, the rest frame duration of the radiation. r= 100, Ep= 1.75 
keV, a= -0.8, l3= -2.4, il=5xl0i3 cm. 
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Figure 9. Symmetric Gaussian intrinsic pulse: pulse width as 
a function of the photon energy for different a, the low energy 
spectral index of the adopted spectrum. Other model parameters 
are the same as in Figure 5. 
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Figure 8. Symmetric Gaussian intrinsic pulse: lag dependence 
on R, the radius of radiation surface. r= 100, E= 1.75 keV, a= 



-0.8, 13= -2.4, t^= 40 s. 



We did not find any peak lag for a rectangular intrinsic 
emission profile, from which a straight-forward conclusion 
regarding the radiation process of the pulse in GRBs can 
be obtained — the radiation intensity must have a decay- 
ing phase in order to produce the observed peak lags. The 
intrinsic decaying phase may be due to the variations asso- 
ciated with hydrodynamic processes, such as the decaying 
of emission caused by density or magnetic field attenuation 
as the shock moves through the shell. 

We have investigated the possible dependences of the 
pulse peak lag upon other physical parameters of the kine- 
matic model. We found: the lag is proportional to the inverse 
of the Lorentz factor; the lag is proportional to the duration 



of the intrinsic radiation t'^; the lag is weakly dependent on 
R; the lag is larger when larger amount of energy is concen- 
trated at Ep (larger a or smaller f}). 

The pulse width decreases with energy {W oc 
^_o.i — 0.2^1^ jjp^ ^ ^ ^j^g observed {W oc E'^-^ ), 

though we found a faster decrease with a larger low-energy 

spectral index a. There must be other energy-dependent 
narrowing mechanisms underlying. Similar to this conclu- 
sion, Dermer (2004) pointed out that other processes in- 
cluding adiabatic and radiative cooling, a non-uniform jet, 
or the external shock process, rather than the curvature ef- 
fect, should be needed to explain the relationship between 
the measured peak photon energy Ep and the measured vFi, 
flux at Ep in the decaying phase of a GRB pulse. They found 
the curvature relationship does not agree with the observar- 
tion (Borgonovo & Ryde 2001). 

For simplicity we have assumed a spherically symmetric 
radiation surface in this paper, though there is some obser- 
vational evidence indicating that the GRB outflow may be 
collimated. Derived from the afterglow observations, the jet 
coming from the GRB central source generally has a half 
opernng angle 6j > F^^ (Frail et al. 2001). Assuming a jet 
geometry with the jet opening angle of 1° or 4° and the 
same parameter values (F, R and td) used in the spheri- 
cal geometry, we calculated the observed pulse shapes and 
the lags; they show no difference from those of the spheri- 
cal geometry. The reason for this is as follows. Even in the 
case of the isotropic radiation surface, the contribution of 
the flux from the outer side region where the observing an- 
gle d is larger than F~^ is relatively very small, because the 
local flux contribution from the radiation surface (i.e., the 
first-part integrand of Equation 5 in Section 2) will decrease 
drastically with as oc (1 + F^e^)-^ when 9 is smafl. 



© 0000 RAS, MNRAS 000, 000-000 



8 R.-F. Shen, L.-M. Song & Z. Li 



ACKNOWLEDGMENTS 

We thank the reforoo for their careful reading of the 
manuscript and for many valuable suggestions. This work 
was supported in part by the Special Funds for Major State 
Basic Science Research Projects of the Ministry of Science 
and Technology, and by the National Natural Science Foun- 
dation of China through grant 10473010. 



REFERENCES 

Band D., Matteson J., Ford L. et al. 1993, ApJ, 413, 281 

Bhat P. N., Fishman G. J., Meegan, C. A., et al. 1994, ApJ, 426, 

604 

Borgonovo L. & Ryde F. 2001, ApJ, 548, 770 

Chen L., Lou Y.-Q., Wu M., et al. 2005, ApJ, 619, 983 

Cheng L. X., Ma Y. Q., Cheng K. S., Lu T. & Zhou Y. Y. 1995, 

A&A, 300, 746 
Dermer C. D. 2004, ApJ, 614, 284 

Fenimore E. E., in t Zand J. J. M., Norris J. P., et al. 1995, ApJ, 
448, LlOl 

Frail D. A., Kulkarni S. R., Sari R. ct al. 2001, ApJ, 562, L55 
Hakkila J. & Giblin T. W. 2004, ApJ, 610, 361 
loka K. & Nakamura T. 2001, ApJ, 554, L163 
Kobayashi S., Piran T. & Sari R. 1997, ApJ, 490, 92 
Koccvski D. & Liang E. 2003, ApJ, 594, 385 
Lithwick Y. & Sari R. 2001, ApJ, 555, 540 

Norris, J. P., Share, G. H., Messina, D. C. et al. 1986, ApJ, 301, 
213 

Norris J. P., Nemiroff R. J., Bonnell J. T. et al. 1996, ApJ, 459, 
393 

Norris J. P., Marani G. F. & Bonnell J. T. 2000, ApJ, 534, 248 
Norris J. P., Scarglc, J. D.& Bonnell, J. T., in Ritz S., Gchrcls N., 
Shrader C. R., eds, AIP Conf. Proc.Vol. 587, Gamma 2001: 
Gamma-Ray Astrophysics. American Institute of Physics, 
New York, p. 176 (astro-ph/0105052) 
Norris J. P. 2002, ApJ, 579, 386 

Norris J. P., Bonnell J. T., Kazanas D. et al. 2005, ApJ, 627, 324 
Piran T. 1999, Phys. Rep., 314, 575 
Piran T. 2004, Rev. Mod. Phys., 76, 1143 

Preece R. D., Briggs M.S., Mallozzi R. S., et al. 2000, ApJS, 126, 
19 

Qin Y.-P. 2002, A&A, 396, 705 

Qin Y.-P., Zhang Z.-B., Zhang F.-W., et al. 2004, ApJ, 617, 439 
Ryde F. & Petrosian V. 2002, ApJ, 578, 290 
Salmonson .J. D. 2000, ApJ, 544, L115 
Sari R. & Piran T. 1997, ApJ, 485, 270 
Schaefer B. 2004, ApJ, 602, 306 

Stern B. in: Poutanen J, Svensson R, eds. High Energy Processes 
in Accreting Black Holes, ASP Conference Series 161. San 

Francisco: ASP, 1999, p. 277 
Wu B. & Fenimore E. 2000, ApJ, 535, L29 



